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Abstract 

We study the localization properties of a disordered tight-binding Hamilto- 
nian on a generic bipartite lattice close to the band center. By means of a 
fermionic replica trick method, we derive the effective non-linear cr-model de- 
scribing the diffusive modes, which we analyse by using the Wilson-Polyakov 
renormalization group. In addition to the standard parameters which define 
the non-linear cr-model, namely the conductance and the external frequency, 
a new parameter enters, which may be related to the fluctuations of the stag- 
gered density of states. We find that, when both the regular hopping and the 
disorder only couple one sublattice to the other, the quantum corrections to 
the Kubo conductivity vanish at the band center, thus implying the existence 
of delocalized states. In two dimensions, the RG equations predict that the 
conductance flows to a finite value, while both the density of states and the 
staggered density of states fluctuations diverge. In three dimensions, we find 
that, sufficiently close to the band center, all states are extended, indepen- 
dently of the disorder strength. We also discuss the role of various symmetry 
breaking terms, as a regular hopping between same sublattices, or an on-site 
disorder. 



I. INTRODUCTION 



An interesting and still debated issue in the physics of the Anderson's localization con- 
cerns the existence of delocalized states in dimensions d < 2, the conditions under which 
they appear, and their properties. This problem, which is, for instance, of relevance in the 
theory of the integer quantum Hall effect |T| , got recently a renewed interest after evidences 
of a metal-insulator transition in two dimensions have been discovered [Q. 

One of the cases in which localization does not occur in any dimension is at the band 
center energy of a tight binding model on a bipartite lattice, when both the regular hopping 
and the disorder only couple one sublattice to the other, i.e. in the so called two sublat- 
tice model Although this is not a common physical situation, its consequences are 

surprising, and seem to escape any quasi-classical interpretation, which, on the contrary, 
provides simple physical explanations of other delocalization mechanisms 0. Already in 
1976, Theodorou and Cohen 0] realized that a one dimensional tight binding model with 
nearest neighbor random hopping has a single delocalized state at the band center (see also 
Ref. 0). Afterwards, Wegner and Oppermann and Wegner showed that a delocahzed 
state indeed exists under the above conditions in any dimension, within a large-n expan- 
sion, being n the number of orbitals per site. Later on, Wegner and Gade ^ argued that 
these models correspond to a particular class of non-linear cr-models for matrices in the zero 
replica limit. They were able to show that the quantum corrections to the /^-function which 
controls the scaling behavior of the conductance vanish at the band center at all orders in the 
disorder strength, thus implying a metallic behavior at this value of the chemical potential. 
Moreover, they showed that, contrary to the standard case, the /5-function of the density of 
states is finite. These results were based upon the non-linear cr-model derived by Gade 
by means of a boson-replica trick method, in a particular two sublattice Hamiltonian with 
broken time reversal invariance. 

More recently, the model without time reversal invariance has got a renewed interest 
for its implications in different physical contexts, for instance models with non-Hermitean 
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stochastic operators, or random Hux models in two dimensions (see e.g. Refs. |1T0|-|T2|). 

In this paper, we present an analysis of a generic disordered tight-binding Hamiltonian 
on a generic bipartite lattice. The starting model is therefore of quite general validity, also 
describing systems with time reversal invariance, and reduces in particular cases to the 
model discussed by Gade or, in the honeycomb lattice, to models of Dirac fermions []T2 



or, finally, to random flux models |IO|,|TT|. By means of a fermionic-replica trick method, we 
derive the generic non-linear cr-model describing the diffusive modes, which we analyse by 
the Renormalization Group (RG). Needless to say, the effective model belongs to the same 
class of non-linear cr-models identified by Wegner and Gade, demonstrating once more the 



universality of this description in the theory of Anderson localization [|T3|. 

Since the work is quite technical, we prefer to give in the following section a short 
summary of the main results. 



A. Summary of the main results 

In this section we shortly present the main results, with particular emphasis to the 
connections and differences with the standard theory of the Anderson's localization. 

We consider a generic bipartite lattice and work with a unit cell which contains two 
sites from opposite sublattices. The Pauli matrices cr's act on the two components of the 
wavefunction, corresponding to the two sites within each unit cell. In this lattice, we study 
a disordered tight-binding Hamiltonian which has the peculiar property of involving only 
both Pauli matrices a\ and a^- In other words, H satisfies the conditions 

{/7,a3} = 0, {i/,ai}^0, {H,o<,} + ^, (1) 

where {...,...} indicates the anticommutator. By means of a path- integral approach within 
a fermionic replica trick method, we find that the low-energy diffusive modes at the band 
center, ii^ = 0, can be represented by the non-linear a-model 

VU{Ry^ ■ VU{Rf 
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S[U] 



dR Tr 



J dR {Tr[uiRrVUiRr]}\ (2) 

where U{R) is unitary and belongs to the coset space U(4m)/Sp(2m), being m the number 
of rephcas. At finite energy E 0, the symmetry of U{R) gets reduced to Sp(2m)/Sp(m)x 
Sp(m), as in the standard case [|1^. The enlarged symmetry is accompanied by new diffusive 
modes which appear in the retarded-retarded and advanced-advanced channels, which are 
instead massive in the standard case. In (Q), axx is the Kubo conductivity (in units of e^/^) 
in the Drude approximation. We find that the new coupling constant, 11, is proportional to 
the fiuctuations of the staggered density of states, i.e. to the following correlation function 

1 ^ e--^^^-^') {{p,{E,R)){p,{E,R')), (3) 

^ RR' 

where the bar indicates the impurity average, and Psi^E, R) is the staggered density of states 
at energy E, 

Ps{E, R)=Y1 MRV^3MR)SiE - e„), 

n 

being (j)n{R) the two-component eigenfunction of energy e„. 

The structure of the above action was derived by Gade and Wegner ^ for a partic- 
ular Hamiltonian. Here, we derive it for a generic bipartite lattice and random hopping. 
Moreover, we provide a simple physical interpretation of 11. 

Going back to (0), Gade and Wegner |^ gave a beautiful proof, based just on symmetry 
considerations, that the quantum corrections to the /3-function of a^x vanish in the zero 
replica limit. In Appendix 0, we outline how their proof works in our case, which is sligthly, 
but not qualitatively, different from the U(N)/SO(N) case they have considered. Essentially, 
one can show that the action (0) posseses an invariant coupling axx + ttlII, which, in the 
m —>■ limit, implies that axx is not renormalized, apart from its bare dimensions. 

On the contrary, both the density of states and 11 have non vanishing /5-functions. In 
d = 1, the system flows to strong coupling, hence we can not access the asymptotic infrared 
behavior. Nevertheless, the starting flow of the running variables indicates that the density 
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of states diverges. In d = 2, 3, the system flows to weak coupling, hence we can safely assume 
that the infrared behavior is captured by the RG equations. Indeed, in two dimensions, the 
density of states diverges a.t E = 0, while, in c? = 3, it saturates to a finite value, although 
exponentially increased in 1/axx- Moreover, 11 has an anomalous behavior in = 2, where it 
is predicted to diverge logarithmically. We explicitly estimate how these quantities behave 
as — i> 0, by means of a two-cutoff scaling approach, as discussed by Gade P]. 

We have also analysed various symmetry breaking terms. The simplest ones are those 
which spoil the particular symmetry Eq.(|ll) of the model at E = 0, i.e. an on-site disorder 
or a same-sublattice regular hopping. These perturbations bring the symmetry of U{R) 
down to Sp(2m)/Sp(m)x Sp(m), as in the standard localization problem. However, by 
evaluating the anomalous dimensions of these terms, we can estimate the cross-over lengths 
above which the symmetry reduction is effective. While in = 1, 2 these terms always lead 
to a localized behavior also at the band center, in c? = 3 the vicinity to the band center leads 
to an increase of the window in which delocalized states exist. 

Finally, if the impurity potential breaks time-reversal symmetry, the matrix field U{R) 
is shown to belong to the coset space U(2m), which indeed agrees with the analysis of Gade 

i- 

The paper is organized as follows. In section ||, we introduce the Hamiltonian. In section 



III, we derive the path-integral representation of the model, by using Grassmann variables 



within the replica trick method, and, in section IV, we study the symmetry properties of the 



action. In section 0, we evaluate the saddle point of the action, while in sections 0, |VII 



VIII| we derive the effective non-linear cr-model describing the long-wavelength fluctuations 
around the saddle point. The Renormalization Group analysis is presented in section 
and the behavior in the presence of on-site disorder, of a same-sublattice regular hopping 
or with broken time reversal invariance is studied in sections and pCIl| , respectively. 

Finally, section pCIV| is devoted to a discussion of the results. We have also included several 



appendices containing more technical parts. 
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II. THE MODEL 



We consider a tight binding Hamiltonian on a bipartite lattice, of the form 

H =Y. Yl ^RR' (cki?' + 4'%) , (4) 
i?eA R'eB 

where A and B label the two sublattices and the hopping matrix elements hji^i are randomly 
distributed. We take a unit cell which includes two sites from different sublattices. In some 
cases, like the honeycomb lattice, this is indeed the primitive unit cell. In other cases, like 
the square lattice, it is not. 

In this representation, the Hamiltonian can be written as 

H = Y. {cl^c^R' + H.c) (5) 

RR' 

where 1 and 2 label now the two sites in the unit cell, while R and E! refer to the unit cells, 
and /i)|fj/ = /i/jz/j. By introducing the two component operators 



Cr 



\ (^2R I 



we can also write the Hamiltonian as 



H = J2^R Hrr' Cr, 

R,R' 

= Y.l {^RR' + ^RR') Cr^iCr' + l {h]^R' - ^RR') 4a2C^,, (6) 
R,R' ^ ^ 

where the cTj's {i = 1, 2, 3) are the Pauli matrices. We notice that, quite generally, the 
Hamiltonian involves both ai and a2, but neither nor ctq, so that it satisfies the conditions 
in Eq. (|1D. 

We can write h^^j^, = t)^^, + r^^, , where t^^p^, are the average values, which represent the 
regular (translationally invariant) hopping matrix elements, while r^^/ are random variables 
with zero average, which we assume to be gaussian distributed with width 



The dimensionless parameter m is a measure of the disorder strength in units of the regular 
hopping. In this way, the Hamiltonian is written as the sum of a regular part, H^^\ plus a 
disordered part, Himp- 

For the regular hopping, we define 



tRR' - - {t^RRl + t^RR^ , WrR' - - {t]^^, - t^j^f^, 



so that the non disorderd part, H^^\ of the Hamiltonian is 

i7(o) = c^R Hfp., Cr, = 4 itRR'(^i + iwRR'(J2) Cr,. (7) 

R,R' R,R' 

Since, for any lattice vector Rq, Irri = t/j+HoR'+iJo' as well as wrr' = wr+r^r'+r^, and, 
moreover, tRRi = Irir while wrri = —wr'r, the Fourier transforms satisfy tk = t\ (tk real) 
and Wk = —wl {wk imaginary). In momentum space, the Hamiltonian matrix, Hj^^ = 
tkO'i + iWkO'2, is diagonalized by the unitary transformation C;, = Ukdj^, with 

Uk = e-'^'^'e''-^"', (8) 

where 

tan^. = ^ = -|^. (9) 

tk Tletf ^ ' 

Indeed, Ulnf^'Uk = tka-s, where el = tl - wl = |4p + 

A. Current operator 

The commutator of the density c^c^ with the non disordered Hamiltonian (|^) is 

^R i'^RR'^l + iWRRl(J2) C^, - C^R, {IrirCJi + iWRRl(J2) 

R' 

or, in Fourier space, 

Y (tk+q - tk) 4f^lCfc+g + ^ {Wk+q " Wk) claiC^+q- 
k 

Therefore, in the long wavelength limit, the current operator in the absence of disorder is 



k 

= H ^fc^fc 4^fe ■ ^Cfc + tkV9kc\B^^k ■ o'Cfc (10) 

A: 

= ^ Vefc dlaadi^^g + CkVOk 4^2C?fc+g, (11) 

the last being the expression in the basis which diagonahzes the Hamiltonian Hq. In (|T0|), 
a = (cTi, (72, CTs), and the vectors 

5fc = (cos6lfe,sin6'fe,0), 5_L,fe = (- sin6'fc, cos6'fc, 0), (12) 

describe intra and inter-band contributions to the current vertex. Notice that the regular 
hopping Hamiltonian can be simply written as H^^^ = J2k ^kc\ ^k ■ <?c^- 

Moreover, since also the impurity part of the whole Hamiltonian, (^, does not com- 
mute with the density operator, in the disordered model the current operator acquires an 
additional term proportional to the random hopping matrix elements, which we discuss later. 



III. PATH INTEGRAL 

The starting point of our analysis is a path-integral representation of the generating 
functional, in terms of Grassmann variables, following the work by Efetov, Larkin and 
Khmelnitskii [|T^. To this end, we introduce, for each unit cell R, the Grassmann variables 
CR]a,p,a and their complex conjugates CR-a,p,a, where a = 1,2 is the sublattice index, p = ± 
is the index of the advanced (+) and retarded (-) components, and the index a runs over 
m identical copies of the model, as in the usual replica trick method. In what follows, by 
convention, the Pauli matrices a's act in the two sublattice space, the r's in the space of 
the Grassmann fields c and c, and the s's in the ± space. 

In order to treat on equal footing both the particle-hole and the particle-particle diffusive 
modes {diffusons and cooperons, respectively), as implied by time-reversal invariance, it is 
convenient to introduce the Nambu spinors \E'ij and defined through 
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Cr 



where cr and cr are column vectors with components CR-a,p,a and CR-a,p,a, respectively, and 

where c = —2x2 is the charge conjugation operator. The action in terms of the spinors is 
[see Eq.(l)] 



S = -J2'^r(e + i^Sg - Hrr,] ^r, 
RR' ^ ^ 



= -Y,^R[E^i^S^- Hf^, )^R> + Y1 Himp,RR' ^iJ' 
RR' ^ ^ RR! 

= 5*0 + Sirapi (13) 

where E are the complex energies of the advanced/retarded components. 



A. Disorder average 

Before taking the disorder average, we notice that, in the spinor notation, 

ArC2R' + H.C. 2^iR^)2R' = 2^2/?' 

Therefore, the impurity part of the action can be written as 

Simp = X! 2r^R/^lR^2R'- 

R,R' 

The generating functional, within the replica method, is 

Z = J ^'^I?^I?rP[r]e"^o-^™^ (14) 

where P[t] is the gaussian probability distribution of the random bonds t^^,. The average 
over disorder changes the impurity action into 



R,R' 

We define 



R,R' 

j2 2w' {t^R^R^y (^iij^2i?') (^2ij'^i/?) . (15) 



Wrr, = 2u' (42^,)' G 7^e, 



so that W* = W-q, and introduce 

where a is a multilabel for Nambu, advanced/retarded and rephca components, and analo- 
gously X2^, as well as their Fourier transforms. By these definitions, 

S.mp = ^ET.W-,Xf^Xt.,. (16) 

This form, as compared to (0), has the advantage to allow a simple Hubbard-Stratonovich 
transformation. Notice that the use of Nambu spinors has the great advantage to involve 
just a single Fourier component of Wrr' . If we write 

where cUg > and (pg = — 0-q, and define 

Simp takes the simple form 

Simp = yYl ^q^hq^i-q = y X! '2'^^ [^1,9^2,-9 + ^2,9^1,-9] 

= ^ E l'^o,qYo,-, - Y,,gYs,.g] , (17) 

^ q 

where we have introduced Yq = Yi + Y2 as well as Y3 = Yi — 1^2- Moreover, our choice of the 
impurity potential, which does not break, on average, the spatial symmetries of the lattice, 
implies that (pq = —9g, see Eq.(|). 

We notice that, if a term is written as 

J2 X-^X^" = X—Tr (X2) , 

a/3 

where X = X"!" and A = ±1, one can always decouple it, by introducing an hermitean matrix 
Q, by the following Hubbard-Stratonovich transformation 
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exp 




= N JvQ exp 


-A 



(18) 

where the normalization factor N~'^ = J VQ exp [— y4^^Tr(Q^)]. In the specific example, 
(17) can be transformed into 



imp 



V q UJq ^ q 

T7 E — [QogQo^q + QsgQs-q] 
^ q ^q 



V 



y^JTr 



Qoq i^cos^Xlq + ^sin^X! 
Qsq fcos^X*_, + ^sin^X*_ 



2 

q vt 



(19) 



If we define Qg = QogCro + iQzqO'z, we obtain 



imp 



9 ^ p,q 



— > 1 r 

V^2UJg 



R ^ p,q 



-q^p+q, 



(20) 
(21) 



where the last term vanishes at g = 0. Notice that the two-sublattice symmetry properties of 
the Hamiltonian, see Eq. (|T]), are reflected in the particular form of the Q- matrix, which only 
contains the Pauh matrices (Tq and a^. In particular, since the tensor Q{R) ~ '${R)'${R), 
(To selects the uniform component while the staggered component of the product of the 
two Grassmann fields. Notice that the electron-Q coupling in the basis can be written 
as a local coupling but for the last term in (pT]) . To simplify the notation of this term we 
find it useful to define the operator L through 



LQ{R) = ^11 e^'"" (l - e^"^'"') Qq- 



(22) 



If we now transform the spinors in Eq.(pOD to the diagonal basis, the coupling to the 
Qq matrix transforms as Up^qQqe^^'''^'^Up. The simplest consequence is that, in the diagonal 
basis, QkiM = Qomm'^o + iQiMM'^T- depends on two wavevectors. However, in the case of 
cubic lattices, see Eq.([AlD, 
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since 0g = — ^^g. Therefore, for cubic lattices, we can also write 

-tY,^RQiR)<^R, (23) 

' q ^^1 R 

where now is the Grassmann field of the dj^ operators, and the matrix Q{R) = Qo{R)ao + 
iQi{R)ai, with Qi{R) = —Qsi^R) defined above. 

Finally, we notice that, in the case of the honeycomb lattice, at the wavevector q^,, con- 
necting the two Dirac cones of the non-disordered dispersion band, Uq^ = 0. This observation 
will turn useful when discussing the long-wavelength behavior of the model. 



S.. 



imp 



IV. SYMMETRIES 

The action (|13]), aX E = u = 0, i.e. at the band center with zero complex frequency, is 
invariant under a transformation "^^i T "^^i if 

cT'c'Hrr,T = HuR'. 

Since the random matrix elements -ff/jj?' involve both Pauli matrices ai and a2, T has to 
satisfy at the same time cT^ aiT = a\ and cT^ a^T = o^- This implies that 

cT*c* =aiT-Vi (24) 
cTiT-Vi =a2T-V2. (25) 



The condition (|25D can be fulfilled only by a transformation T = TqGq + Tads 



Under such a transformation 

Q cT'c'QT = aiT-^aiQT = asT^VsgT 
Since aiQai = Q\ then 
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Hence, the transformation is also unitary, = T^^. Moreover, such a transformation leaves 
the Q-manifold invariant, which implies that our Hubbard-Stratonovich decoupling scheme, 
which makes use of Q = Qqc^q + iQ^cr-i, is exhaustive. 
The unitary transformation, T, can be written as 



exp 



(26) 



where 



Wl = -Wo, = -Ws. 



im- 



cwy 

cW^c' 



In addition, we must impose the charge conjugacy invariance, which, through Eq. 
plies that 

-Wo = wl 

W3 = -wl 

The number of independent parameters turns out to be 16m^, which suggests that T is 
related to a unitary group, specificaly U(4m), as argued by Gade and Wegner §]. In fact, 
we can alternatively write 



T 



e 2 



v 









e 2 



U 



cWd 



J 



(27) 



where U is indeed a unitary transformation belonging to U(4m). The invariance of (|TB[) 
at finite frequency, u ^ 0, implies the additional condition cT^c^s^T = S3, which reduces 
the number of independent parameters to 8m^ + 2m, lowering the symmetry of U down to 
Sp(2m). 

li E ^ 0, T has to satisfy also 



cT'c'T = oiT-'oiT = 1. 

This implies that, at finite energy, i.e. away from the band center, T does not contain 
anymore a as-component. Indeed E lowers the symmetry of U down to Sp(2m), which is 



further reduced to Sp(m) x Sp(m) by a finite frequency, as in the standard situation |^ 
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V. SADDLE POINT 



The full action 



k,q 

by integrating over the Nambu spinors, transforms into 



(28) 



-Tr In 
2 



E + z^S3 - +iQ~ iLQ 



(29) 



The saddle point equation for homogeneous solutions at £" = and u; = O"*" is 

Q = ^^J^ (^0+S3 - efc + iQ) + [lO+ss + ek + iQ) , 
where uq = J2r-r' '^u^iiR-wY ■ The general solution is 

TT 

Qsp = ^^^0P(0)S3 = SS3, 



(30) 



with p(0) being the density of states at = 0. In order to distinguish transverse from 
longitudinal modes, it is convenient to parametrize the Q-matrix in the following way 

Q{R)p = aiT(/?)-Vi [Qsp + P(i?)] TiR) = Q{R) + a{r{R)-^<jiP{R)T{R). (31) 



Here P{R) describes the longitudinal modes, which we discuss more in detail in section |VT|, 
and T the transverse modes. Namely, T has the form given in Eq. 



T{R) = exp 



exp 



— — ^0 H — 0-3 



(32) 



with exp[(iyo + ^3)72] belonging now to the coset space U(4m)/Sp(2m). This amounts to 
impose that 



{1^0,5} = 0, [W3,s]=0, 



by which it derives that 
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/ 









(.A 
















V 


-St 








iC ^ 



Q{R) = aiT{R)-'aiQspT{R) = g,^e^o(«)-o+^3{ii)a3. 
In the ± space, we can write 

Wo 

where — A, — C, and additionally, since cW^c^ — giW^gi — 
A, cC*c* = C and cS*c* = St. By writing 

A = AoTo + i (Ain + A2T2 + A3T3) , 
B = BoTo + i (Sin + B2T2 + SsTa) , 
C = CoTo + i (Cin + C2T2 + CaTa) , 

we find that the above conditions imply that, for i = 0, . . . , 3, 

Bi,Ai,Ci e Tie, 

and 



(33) 



(34) 

-GiWai, then cA^c^ — 

(35) 
(36) 
(37) 



(38) 



while, for j = 1,2,3, 



^0 — ^0 — C'o) 



— -A*, — -C*. 



(39) 



(40) 



VI. EFFECTIVE ACTION 

In this section, we derive the effective field theory describing the long wavelength trans- 
verse fiuctuations of Q{R) around the saddle point. In the case of honeycomb lattices, we 
should worry about the momentum component of Q which couples the two Dirac cones. 
However, one can see that the free action of Q diverges at this wavevector, so that we are 
allowed to ignore the fiuctuations around this momentum. 
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A. Integration over the Grassmann fields 



As we said, by integrating over the Grassmann variables, we obtain the following 
action of Q: 



H — Tr In 
2 



E + i^s; - +iQ- iLQ 



(41) 



We start by neglecting the longitudinal fluctuations. Then, since Q = T^'QspT, where we 
define T = aiTai = a2Ta2, we can rewrite the second term of S[Q] as 



^Tr In (eTT^ + i^fs^T^ - TH^^'^T^ + iQ,p - , 



where we define 



(42) 



V = iTLQTl 



Since H^j^, involves either ai and (T2, while T involves uo and as, then 

= r(i?)t4°), - vf (i?)t . (r - R') hSI, 

r(0) 



+ -d,,T{Ry {R, - i?:) (R, - R'A HZ- 



Therefore the term TH^^^T"^ which appears in (H^ can be written at long wavelengths as 



f{R)H^^l,T{Ry = hZ - ^f{R)Vf-\R) ■ J '^'\R')nR' 



AO) 



+ lf{R)d.,fiR)-' (R, - R[) {r, - i?;.) 4% 



= ^Pj^, + UrR> 



(43) 



where we used the fact that the long-wavelength part of the current operator in real space 



is, see (p3D, 



JiR') 



R 



Notice that in (^) only the current vertex which derives from the regular hopping appears, 
which is not the full current operator. 
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Moreover, a further current-like coupling will arise from the expansion in V (see below). 
To this end, in the long- wavelength limit, the operator L, see (^2]), can be approximately 
written as 



LQ{R) ^ -/3 ■ VQ{R)a, - i (/3 ■ Q{R), 



(44) 



where = V0(g = 0)/2. 

Having defined U, (H2]) can be written as 



^Tr In (^ETT^ + I'^fssT^ -U-V- H^"^ + zQ^p^ 
= -^Tr In G + ^Tr \n(l + G ETT^ + Gi'^fs-iT^ - G f/ - G , (45) 

where G = {—H^^^ + iQsp)~^ is the Green's function in the absence of transverse fluctuations. 

The effective field theory is then derived by expanding S[Q] up to second order in U and 
V, and first order in E and u. In this way we get the following terms. 



B. Expansion in the Q free action 



The free part of the action 
can be expanded at small q. Since iVg = oj^„, then 



leading to 



1 



2ujQ 
1 



dRTr 
dRTr 



Q{R)Q{R)^ 



+ 



7 



2Vuoo 



Q 



+ 



7 
2ujQ 



dRTr [vg(i?) ■ vg(i?)^ 



(46) 



The second term is a contribution to the current- current correlation function of the part of 
the current vertices proportional to the random hopping. 
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C. Expansion in E 



Expansion of (E^) in E gives 



f Tr (GTT^) = -i—Tr (QspTT^) = -i—TrQ. 



D. Expansion in lo 

Expansion in lo gives 



i-Tr ^GfsT^) = —Tr{ssQ) 



E. Expansion in U 

The second order expansion in U contains the following terms: 

1. 



^TriGU) 



and 



- ^Tr (GUGU), 

Taking in (^91), the component of U containing second derivatives, we get 
~Tr {f{R)d,^f{R)-' {R, - R!,) {Rj - R'^) H^GiR', R)} 
By means of the Ward identity ( |B6D , the above expression turns out to be 

_^Tr{f{R)d,,f{Rr'}, 
which, integrating by part, is also equal to 

y + + 

_ ^j,^ {f{R)dif{R)-^f{R)djf{R)-^} 
- Y^xtj^Tr [DiDj + DiSscriDjS^ai] . 



Here we have introduced a matrix D{R) with the i-th component 



D,{R) = Do,i{R)(ro + Ds,iR)a3 = f{R)dif{R)-\ 



(54) 



The part of (H^) which contains first derivatives gives rise to a boundary term 



dRTr 



VW{R)(T^ 



(55) 



where S has been defined by Eq.(^), which we discard by taking appropriate boundary 
conditions. 

Let us now analyse the term (|50|) , where we have to keep of U only the part containing 
first derivatives. By making use of (^), this term is, in momentum space, 



kg 
1 



TVT 



-1 



mG{k + q) 



kq 



Jk^'^m 



k R 



(56) 



valid for small q. Here the matrix jj^^^ = VktkO'i + iV kWk<^2- The non vanishing terms in 



(0) have both D's either Do^o or -Daas [see (^4])]. 

In the diagonal basis, upon defining, as we did in Eg. (|30|) , Qsp = Sss, with S 
7ru;op(0)/4, the Green's function is 

1 



G'o(/i;)(ToS3 + G3{k)a3So. 



(57) 



Going back to the original basis. 



Gik) = UuQ{k)Ul = Goik)aoS3 + G3ik)B, ■ as^, 



(58) 



where Bf^ has been defined in (|T2|). Therefore, 
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(T3G'(A;)a3 = Go{k)aoSs - G'i{k)Bk ■ aso = -siG{k)si. (59) 
By means of (|59D , we find that 

asJ.^'^Gikyas = -J,^'^ a^Gik^a^ = J,^''^G{k)-, (60) 
from which it derives that (^61) can be written as the sum of two different terms 



^X^j-Tr [D^Dj - D,s^a,DjSsa,] (61) 
+Y^xti'^Tr [D,D, + DiS^aiD,s:,a^] . (62) 



By summing (|1]), (|2|), (H) and (0), we get 

1 



16 



which is equal to 



2tt 



■ajf j dRTr [diQ{R)d,Q{R)^) , (63) 



where 



.(0) _ !_ ( 



is the Kubo conductivity with the current vertices which involve only the regular hopping 
[cf. Eq.(i3. 

F. Expansion in V 

The expansion in V up to second order, gives two terms 

-\Tr{GV), (64) 

and 

-^-Tr{GVGV), (65) 
20 



In addition, we must also consider the mixed term 



- -Tr (GUGV) 



(66) 



The first order term 



IS 



-]-Tr {GV) = —Tr {Q,pV) 
— [dRTr (Q{Ryp ■ VQ{R)a3 



1 



2tu, 



dRTr (p-VQiRYP- VQ(i?)) . 



(67) 
(68) 



The first term ( p7D is another boundary term, which we neglect. The second order term, 
Eq.(|65l), gives 

(69) 



- ^ E {Goik)' + Gsik)') J dRTr [0 ■ VQ{R)^P ■ VQ(i?) 
Notice that the saddle point equation implies that 

V k ^0 
By using the above equation in (|68|) , we find that (|69|) plus ( |68D give 



-— Gsikf j dRTr [p ■ VQiRYP ■ VQ{R) 



1 



dRTr 



P-VQiRyp-VQiR) 



(70) 



The contribution of the mixing term, (pB|), can be evaluated in a similar way, giving 



1 



—y 



dRTr 



Vek-VQ{Ry(3-VQ{R) 



(71) 



We notice that, since 0^ = —9q, then ([70|) and (|7lD can be included in (|63D if the following 
redefinition of the current vertex is assumed: 



=VekBk-a + ek {vOu - VOo) B^^j, ■ a. (72) 

Therefore the Kubo conductivity which appears in (^) has to be calculated with the above 
expression of the regular current vertex. This has notable consequences. First of all, for 
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cubic lattices, the interband contribution vanishes, hence the Kubo conductivity coincides 
with the true one, since the enlargment of the unit cell was artificial. This is compatible 
with (p3D, where we showed that the impurity action is local in the basis which diagonalizes 
the regular hopping, implying that the regular current vertex contains only the intraband 
operator. 

To conclude, the effective action so far derived is therefore 
S[Q] = 7T- f dRTr \Q{R)Q{R)^^ 

27r 



+ / dRTr {d,Q{R)dMRy 

+ I dRi—Tr {Q{R)) - ^Tr {s;Q{R)) . (73) 

J UJq 2UJq 



VII. LONGITUDINAL FLUCTUATIONS 



The full expression of the Q-matrix we must indeed consider is the one given by (pTD, 

Q{R)p = f{R)-' [Qsp + P{R)] TiR) = Q{R) + f{R)-^P{R)T{R) 

^Q{R) + S{R), (74) 

where T[R) involves transverse fluctuations and 



P(R) = (Pooso + P03S3) (To + i (P31S1 + 



(75) 



being all P's hermitean. Charge conjugation implies that cP*c* = P. The fleld P{R) takes 
into account longitudinal fluctuations which are massive. Let us deflne r[R — R') the Fourier 
transform of uj^^. Then, the free action of the Qp fleld is 
1 



S[Qp] = ^ I dRdR'T{R- R')Tr Qp{R)Qp{R!)^ 



2a;n 



dRTr 



Qp{R)Qp{Ry 



-- I dRdR'r{R-R')Tr (Qp{R) - Qp{R!) ) (Qp{R)^ - Qp{R!^^ 



(76) 

(77) 
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Since ggt = Q2^^ gives 

/ dRTr \P{R)P{Ry + 2Q,pP{R) + Q\ 
The second term cancels with the first order expansion of Tr In Gp, since Qsp is the saddle 



point solution. What is left, i.e. 

1 



dRTr 



P{R)P{Ry] , (78) 



2uJo 

is actually the mass term of the longitudinal modes, since the second order expansion in P 
of Tr In Gp is zero. The other term, (^), can be analysed within a gradient expansion of 
QpiR) - Qp{R!) = VQp{R) ■ (r — R'^ + . . .. The details of the calculations are given in 
Appendix P 1| , so that, in this section, we just present the final results. 
The free action of the longitudinal fields is found to be 

S.lP]-^E^TriP,Pl). (79) 

Here, we neglect the the contribution of the invariant measure, which, in the zero replica 
limit, gives rise to fluctuations smaller by a factor than Eg. ([79D . The integration over P 
with the above action has several important consequences for the action of the transverse 



modes (see Appendix |C1 



First of all, all the terms of the Kubo conductivity with the random current vertices are 
recovered. In addition, we find a new operator 

27r 



U J dRTr [Q^{R)VQiR)a3] ■ Tr [Q^ {R)VQ{R)a3\ , (80) 



8 ■ 32S4 

which has contributions from two different terms. The first one is obtained by expanding 



each Green's function in ([50| ) at first order, Gp = Gq — iGqPGq, and the second is derived 
by ([771). They are analogous to the components of the Kubo conductivity with regular and 
with random current vertices, respectively. 
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VIII. EFFECTIVE NON LINEAR cj-MODEL 



In conclusion the final expression of the action of the transverse modes in the long- 
wavelength limit is 

2n 



S[Q] = ^a,.,. j dRTr (VQ{R) ■ VQ(i?)t) 
+ / dRt—Tr (QiR)) - ^Tr {s,Q{R)) 

- ^:|^n J dRTr [Q\R)VQ{R)as\ ■ Tr [Q^ {R)VQ{R)as\ , (81) 



where we make use of the fact that, in the models we consider, atj = Sijaxx- Since Q{R) = 
QspT{RY, see Eq.(^), expressing T{R) by means of U{R) as in Eg. (pTf) , the action at 
E = uj = can also be written as 

27r(T,.T 



S[U] = / dR Tr \VU{R)-^ ■ VU{Rf 

16 J L 

I dR {Tr[uiRrVUiRr]}\ (82) 



as anticipated in the section [1 A| . As compared to the non-linear cr-model which is obtained 
in the absence of the particle- hole symmetry [Q, the above action differs first of all because 
of the symmetry properties of the matrix field U{R), which describes now the Goldstone 
modes within the coset space U(4m)/Sp(2m). Moreover, it also differs for the last term 
of (0), which, in the general case, even if present, is not related to massless modes. An 
analogous term was originally obtained by Gade P] in a two sublattice model described by 
two on-site levels with a regular hopping of the form ifijj?' = ti?_R'Cri, and a local time-reversal 
symmetry breaking random potential Himp,R = Wircti + W2r<J2, which we discuss in section 
pOT , 



Although the action may be parametrized by a simple unitary field U{R) as in (|8^), we 
prefer to work with the matrix Q{R) which has the more transparent physical interpretation 
Q{R) ~ ^{R)'^{R). 

Finally, it is important to notice that either axx and 11 have contributions from both 
the regular and the random current vertices. This implies that, even in the limit of strong 
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disorder, in which the average hopping is neghgible with respect to its fluctuations, these 
constants are flnite and become of order unity [0,|lT| . 



A. Gaussian Propagators 



At second order in W, the dispersion term 



32S2 
32 



dR TriVQ'^VQ 



27ra, 



dR Tr (VWVW) 
J dR ATr (VBVB^) + 2Tr (VAVA + VCVC) , 



where A, B and C are deflned through Eqs.(p5|), (|36D and (|37D . For the -B's we flnd the 
quadratic action 



so that 



1=0 k ah 



{BiAk)BjA-k)) = 6ij6acSbdD{k), 



(83) 



where 



D{k) 



1 1 



For the A's we have to take into account also the disconnected term: 



27rn 



32 ■ 8S4 
27rn 



dRTr 



Q\R)VQ{R)a3\ ■ Tr \Q\R)VQ{R)cr3 



64 
27rn 



dRTr 



■ Tr 



dRTr 



dRTr 



VA + VC 
VAo + VCo 



■ Tr 

Tr 



VA + VC 
VAo + VCo 



(84) 



The non vanishing propagators are 
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(Ao,afe(fc)Ao,cd(-fc)) = D{k) {SacSbd + ^adhc) 



D{k)- 



n 



ab'Jcdi 



an,x + nm 

(Co,ab(fc)C'o,cd(-fc)) = D{k) {6achd + ^adhc) 



D{k) 



n 



{Ao,abik)Co,cdi-k)) = -D{k) 



n 



ab'Jcdi 



ab'Jcdi 



where m is the number of rephcas, while for i = 1, 2, 3 

{Ai^ab{k)Ai^^d{-k)) = D{k) {5ac5bd - ^adhc) 
{Ci^ab{k)Ci^cd{-k)) = D{k) {5ac5bd - ^adhc) 



(85) 

(86) 
(87) 



(88) 
(89) 



Notice that the particular symmetry of the two-sublattice model leads to additional 
diffusive modes in the retarded-retarded and advanced-advanced channels, which are not 



massless in the standard case 16 . 



B. Physical Meaning of IT 

Let us introduce an external source which couples to the staggered density of states, 
which is accomplished by adding to the action a term 

j dR'^RS5a3XiR)^R, 

where, in the replica space, the source A^^^ = XaSa,i3- The fluctuations of the staggered 
density of states is obtained by the derivative of the partition function with respect, for 
instance, to Aq, and Xp, with a ^ (3. Inserting the source term in the action, and integrating 
over the Grassmann fields after introducing the matrix Q, leads to the following expression 
of the staggered density of states fluctuation in terms of Q: 

F{R, R') = -^{Tr [QaaiR)s3cr3] Tr [Q^f,iR')s^sa3]). (90) 

The gaussian estimate of the above correlation function at momentum k is given by 
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F{k) = -IQ {[Ao,aaik) + Co,aaik)] [A^M-k] + CoM-k)]) 

7T LOq 



uJq an,x + nm 



(91) 

Therefore, 11 is directly related to the singular behavior of the staggered density of states 
fluctuations. 

IX. RENORMALIZATION GROUP 



In this section, we will apply the Wilson-Polyakov Renormalization Group (RG) proce- 



dure 1 14, 16 1 to analyse the scaling behavior of the action. Indeed, some of the calculations 
which we present are redumdant, given the proof by Gade and Wegner that the /3-function 
is zero |^ (see Appendix]^). Nevertheless, other results besides the conductance /3-function 
are important, so that we describe the whole RG procedure. 



A. RG equations 



In the spirit of Wilson-Polyakov RG approach [0,0 , we assume that 

T{R)=Tf{R)Ts{R), 

where Tf involves fast modes with momentum q G [A/s,A], while Tg involves slow modes 
with q G [0,A/s], being A the higher momentum cut-off, and the rescaling factor s > 1. 
Within an e-expansion, where e = (i — 2, we define 



A dk 



:D{k) = L 



Ia/s {27rY 47rV, 
It is straightforward to show the following result 



Ins + 0(e). 



Tr 



Tr 



VQ\-VQf 



+2Tr 



Ds(TiQfDsQ\ai 



2T?Tr 



DM. 



+4Tr DMVQf 



(92) 
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where Qf = TlQgpTf and Ds = TgVTj'. Moreover, 



-Tr 



Tr 



Q^VQa; -Tr Q^VQas 



VWs + VWf) as] ■ Tr [(VW, + VWf) 



(93) 



Since the fast and slow modes hve in disconnected regions of momentum space, only the 
stiffness term (p^ ) generates corrections. By expanding the terms coupling slow and fast 
modes up to second order in Wf, the stiffness generates an action term for the slow modes 
which, after averaging over the fast ones, is 



where 



27rCTa;; 

32S2 



dRTr 



1 



(94) 



271 a ^ 



dR - 2Tr 



32S2 

+2Tr [DaiQspDQspWfa, 



Da,QspWfDWfQsp(T, 



(95) 



and 



^ 27ra^ ^^^^ 
32 



DWfVWf 



(96) 



The explicit evaluation of these terms is outlined in Appendix 0. Here we just give the final 
results. The Kubo conductivity is renormalized according to 



CTxx (1 -4Lm)a^ 



(97) 



while the 11 factor 



n n + 4La^ 



(98) 



For what it regards the renormalization of E and u, we notice that 



Q — Q spTsTj-Tg 



QspT^ + QspTsWfTs + -QspTsW]Ts. 
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Since the slow and fast degrees of freedom are defined in different regions of momentum 
space, only the second term is relevant. By means of Eq.( pi5| ), we find that 

«^«.{l + i(2-8™ + ^-!L_)L} (99) 

This leads to similar corrections to E and u, which will have the same scaling behavior. 

Finally, to describe the cross-over behavior in the presence of symmetry breaking terms, 
we also need the scaling behavior of the operator Tr [Q{Ry]. We find that 



{TrQ^)f = {Tr Q spT^TfTsQ spT^TfT, ) 



1 + 2 ( 2 - 4m + ^ 



TrQ'i - L {TrQ.f . (100) 



a^^ + mil. 

To implement the RG, we have to rescale the momenta in order to recover the original 
range [0, A]. This is accomplished by the transformation k ^ k/s, or, equivalently, R Rs. 
Therefore, the stiffness as well as the fluctuation terms acquire a scaling factor ~ 1 + elns, 
while the E and u terms a factor s"^. Hence, after defining t = l/(47r^(Txx), c = l/(47r^n), 
and A the coupling constant of the operator TrQ"^, we get the following /3-functions 

^^ = -et + 4mt2, (101) 

= -ec - 4c2, (102) 

& = d£ + 4(2 + ^), (103) 

(3x = dX + 2Xt (2 + — - — ) . (104) 



c + mt 



(105) 



At finite energy, E ^ 0, we may use a two-cutoff scaling approach 0. Namely, we can 
follow the previous RG equations up to a cross-over scale, Scross = s{E), at which the energy 
has flowed to a value Eq of order S, which plays the role of the high-energy cut-off in the 
theory. Above this scale, we must neglect all contributions coming from the W3 modes, 
which acquire a mass. That is, we must abandon the RG equations ( p.01| )-( p!0^ , and let the 
coupling constants flow in accordance with the standard RG equations, which amount only 
to a renormalization of t according to 
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A = -et + t\ (106) 

If, by integrating ( [L06| ), the inverse conductance t{s) flows to infinity, signalling an insulating 
behavior, then we can define a localization length ^ioc{E) as the scale at which t has grown 
to a value of order unity. 



B. RG in d = 2 

In d = 2, the solution of the RG equations for t, c, and E is 

t{s) = t{l) = to, 

= ^TT + 4 In s = — + 4 In s, 
c{s) c(l) Co 



In 



E 



2 V Co 



1 

2' 



In s + -to l'^^ "5 



At finite energy, the crossover length, s{E), to the standard, non particle-hole symmetric, 
model, which, as previously discussed, is defined through E[s{E)] = £^0 ~ S, is given by 

1 



s{E) = exp 



2A 



'b^ + AAln^ - B 



(107) 



being 



Above s{E), t flows according to Eq.( |106| ) with e = 0, hence it grows to infinity, implying 
that the wavefunctions are localized for any E ^ 0. The localization length C,ioc{E) ~ s{E), 
apart from a multiphcative factor which is ~ exp[(l — to) /to] if to <^ 1. We see that, for 



2 + ^f2 + ^ 
2 V Co 



the localization length has a power law behavior, namely 

'Eo\^ 



iioc{,E) oc 



E 



otherwise, at very small energy, it diverges slower. 
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The density of states renormalizes like 

ln^ = ^f2 + ^Vn^ + -^oln^-. 
Po 2 V Co/ 2 

At finite energy, the density of states flows until s < S{E), after which it stays constant. 
This implies that the renormalized value is obtained by 

ln^ = l„^^ = ln^-21n.(B). 

Po Po E 

leading to 

the last equality valid at small energy. 



C. RG in d = 3 



In d = 3, 



t{s) = toS \ 
c{s) 



hence both t and c running variables vanish for s — > oo. The cross-over length diverges, as 
we approach E — 0, approximately like 

(109) 

Hence, the density of states 
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saturates at = to a value exponentially increased in to with respect to the bare po- 



At E ^ 0, the inverse conductance above Scross = siE) flows according to ( |106| ) with 
e = 1 and boundary condition t{scross) = to/scross- We find that, if 

'^O ^ 'Across J 

the system is metallic, otherwise it is insulating, with a localization length 

This results implies that, for any amount of disorder, sufficiently close to = 0, all eigen- 



functions are delocalized, in agreement with recent numerical results However, if the 
disorder is gaussian, as we assumed, the random hopping model with zero regular hopping 
seems to be characterized by an inverse Drude conductance, to, which is an increasing func- 
tion of \E\, being smaller than the critical value e at E = (see also Ref. 0). In this case, 
the presence of a finite mobility edge in ci = 3, even for zero regular hopping, would not 
depend crucially upon the intermediate RG fiow in the vicinity to the band center. Never- 
theless, we expect that to a.t E = varies for different kinds of disorder, and eventually it 
may become greater than unity. In this case, it is just the vicinity to the band center which 
makes it possible a finite mobility edge. 

X. ON SITE DISORDER 

In this and in the following section, we analyse various symmetry breaking terms, which, 
in the two sublattice representation, contain ctq and a^, hence spoiling Eq.(||). 
We start by adding an onsite disorder 

oSimp = I^Ui^RCIrCi^r + U2,rC'2,rC2,R = I^c'r ^ ^^0 + ^ (^3 j Cr, 

R R 

where (Mi,^) = and {ui^RUj^R') = SijSRR'v"^. Within the path integral, this term becomes, 
once average over disorder is performed. 
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<J>Jimp — „y 2^ -< l,g -'1,-5 "I" -'2,q -'2,-g 
9 



V 



^0,q "T -'3,q -'^3 



3,-g- 



By adding this term to (jT^), we get 



1 



Simp + SSimp = JTT ^{^q + v'^)Tr (Yo.g^O -9) 



3-q) 



111) 



If we assume that the onsite disorder is weak, i.e. Uq > v"^ at small q, the consequence is 
that the Q free action becomes 
1 



"! „,2 ^^ [Qo,qQo,-q] + —^—l^Tr [(53,g<53,-gJ 



-Tr 



QqQ\ 



(112) 



Therefore, the on site disorder introduces a mass in the propagators. Specifically, since 
2iQ3(J3 = Q — , the mass term can be written as 



4K - 



Tr 



Q, - Ql) Q, - Ql 



Close to the saddle point, QQ'^ = Q'j.p, and, for small q, we get 



dR Tr 



Q{R)Q{R)+Q{R)'^Q{R)^ 



which, at second order in W , reads 



_2fY?_ 



dR Tr [W{R)W{R) + W{R)s2W{R)s^] 
dR Tr [W^{R)W^{R)] . 



(113) 



In the presence of this term, we could proceed, as before, in the framework of two cutoff 
scaling theory. That is, we apply the previous RG equations until the above term becomes 
of the order Spo, i-e. up to the scale which, by Eq.(|104|), is 



exp 



{2 A 



?2 + In A - 



(114) 
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in (i = 2, where A = 2t'^ and 5 = 2 + to(2 + ^o/co), while Scross — ^^^'^ in c? > 2, where 

A oc — , 

in the hmit of small v. Above this scale, the propagator gets fully massive, and the inverse 
conductivity flows with the RG equation (|106| ). In c/ = 2, this implies that, ultimately, the 
system gets localized, although the density of states has increased in the first stage of the 
RG. 



XI. SAME SUBLATTICE REGULAR HOPPING 

We can also introduce a particle-hole symmetry breaking term, by adding to the Hamil- 
tonian a regular term connecting same sublattices, e.g. 

Expanding the action in 5H, after integrating over the Nambu spinors, we get an additional 
term 

(115) 
We define 

uJo ^ k 

so that (|115|) becomes 

- z— Tr [Q (Aodo + Agas)] . (116) 

The Ao-term acts like an energy term. This implies that, if we just shift the chemical 
potential, we do recover the same scenario as in the absence of this term and at ii^ = 0. On 
the contrary, the As-term is always a relevant perturbation, whose strength increases under 
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6S[Q] 



-Tr 



fSHT^G 



RG iteration as the energy E. We can define a crossover scale s cross, which has the same 
expression as s{E) in ( |107| ) and ( |109D , for = 2 and d = 3, respectively, provided E A3. 
Above this scale, the W3 modes get fully massive and their contribution to the RG flow 
drops out. 

(3) 

Sometimes t]^j^, = 0, as, for instance, for next-nearest neighbor hopping in a square 
lattice. In this case, A3 = and we need to evaluate the second order term 

1 



6S[Q] = -Tr 



f6HT^Gf6HT^G 



If we define F{R) = T{R)T{R)\ which contains either o"o and a^, this term can be written, 
at long wavelengths, as 



By introducing. 



^S[Q] = W{tfyTr[F,G{k)F_,G{k)] 

^ k 



n' = \^Y.{if)'Tr[a,G^a,G'^ 



where p,q = ±, the following results hold 



1 4C 1 4C 

2 ujo 2 uJo 



where C is a constant of dimension energy square, with order of magnitude given by the 
typical value of (^1°'*) close to the surface corresponding to = 0. Therefore we can write. 



SS[Q] 



C 

2uJn 



dRTr 



F{R)FiR) - F{R)s,a^F{R)ssa, 



const. + 



G 



2uJc^T? 



dRTr 



Q{Rf 



(117) 



which is a mass term for the propagators, similar to ( |113|) . Therefore, a same-sublattice 
hopping introduces a cross-over length analogous to ( |114| ), with 

Aoc-. 

Above this scale, the contribution of the modes to the RG flow has to be dropped out. 
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XII. TIME-REVERSAL SYMMETRY BREAKING 



If the random hopping breaks time-reversal symmetry, i.e. 

RR' 

with both real and imaginary part of r^^, gaussian distributed, after averaging, the impurity 
action can be written as 

Simp = 77 X! ^-qTr [Xifi^gX2fi-q + -^l,3,g-^2,3 -g] ; (US) 
Q 

where 

KIr = ^iR^o^m^ KIr = ^IR^A, 
with the indices a and (3 running only over the replicas and the advanced/retarded compo- 
nents. This implies that the manifold in which Q varies contains in this case only tq and 
r3 components. Indeed, as in the time reversal invariant case we are able to parametrize 



the 8m x 8m matrix T in terms of a 4m x 4m matrix U G U(4m)/Sp(2m) [see Eq.(|27|)], 
similarly, without time-reversal symmetry, T can be parametrized by means of a 2m x 2m 
matrix U G U(2m), in agreement with Gade The effective non-linear a-model is not 
modified, but the expressions (|D1| ), ( p2| ), (|D3|), (|D4|) , and (P5| ) have to be substituted by 

{BatPbcBl) = 2D{k)6adTr (Pq) , (119) 

{BabPbcBcd) = 0, (120) 

= 2D{K)5adTr (Po) - D{k) ^^Pad, (121) 

cr^.^ + Ilm 

(a^ncCcd) = 2D{K)5adTr (Pq) - D{k) ^^Pad, (122) 

cr^r^ + Ilm 

{A,,P,cCJ) = -D{k) ^^Pad. (123) 

a^^ + Urn 

where P = Pq + iPsT^. Hence, the RG equations at m = are, in this case, 

A = -et, (124) 
pc = -ec - 2c^ (125) 
PE = dE + ^, (126) 

36 



which coincide with those obtained by Gade 0. 



XIII. DISCUSSION AND COMPARISON WITH THE STANDARD 

LOCALIZATION THEORY 

In this section, we summarize the main differences between the model and the stan- 
dard non-linear cr-model which is derived in the theory of Anderson localization |T3|,ITB 



placing particular emphasis on the properties of the Q-matrix. The specific form of the 
off-diagonal disorder we consider, which only couples one sublattice to the other, leads, via 
the Hubbard-Stratonovich decoupling, to the introduction of a space- varying 8m x 8m com- 
plex Q-matrix, Q = QoO"o + ^Qs^s- Here, Qo and Qs are 4m x 4m hermitean matrices, 
of which matricial structure refers to the retarded/advanced, spinor particle/hole and m 
replica components. Contrary to the standard case, Q is not hermitean. 

The evaluation of the saddle point, Qsp = cqTos^ (section^), as well as the derivation of 
the efective action (section |VI|) are analogous to the standard case [T^|T^. (We recall that 



(Tj, Si and Ti indicate the Pauli matrices, including the unit matrix, acting on sublattice, 
advanced/retarded and spinor components, respectively.) The non-linear a-model, Eq.(pi|), 
is obtained by integrating out the longitudinal massive Q-fluctuations and only keeping the 
transverse soft modes. The real novelty with respect to localization theory is not in the 
structure of the effective action. Indeed, the new term in (|8ll) , namely 



n / dR 



Tr {Q{RyVQ{R)a3 



even if present, would be irrelevant in the standard case. On the contrary, the essential 
difference, as expected, lies in the ensamble spanned by the soft modes at the particle-hole 
symmetry point E = 0. We get Qsoft = T~^QspT, where the unitary matrix T only contains 
(To and (T3, 

T = exp [ 

and 
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T = aiTai = O'lTa^ = exp 



2 

These expressions derive by the conditions (ffl), which fully specify the model, as shown in 



section f^. In that section, we also showed that the ensamble can be expressed in terms of 
unitary 4m x 4m matrices 



U = exp 



2 

as argued by Gade and Wegner g]. Selecting the subset which leaves the saddle point 
invariant gives Eq.(0) with U G U(4m)/Sp(2m). In terms of T, the condition T~^QspT ^ 
Qsp, leads to the requirements [Wq, S3] 7^ and {W^, S3} 7^ 0, which implies that Wq is 
off-diagonal in the energy retarded/advanced space (as in the standard localization theory), 
while Ws is diagonal. In other words, the omogeneous and staggered modes, Wq and W3, 
respectively, have different structure in the energy space. The energy diagonal iy3-modes 
betray the presence, at E = 0, of diffusive poles in the disorder averaged products of retarded 
and advanced Green's functions, GrGr and GaGa, with Gr^a = {—H±iO^)~^. In the 



localization theory [|T3|,|l6|, only the mixed products GrGa have a singular behavior. This 



explains why singular corrections to the density of states (which involves connected diagrams 
with same energy Green's functions) are present in the two-sublattice model, while they are 
absent in the localization theory. 

In a square lattice, the energy diagonal modes have the transparent meaning of density 
fluctations with wave- vector q nearby the nesting vector G = (7r,7r, . . .), see Appendix A. 
Indeed 

R&A 

= ^ e-^('^+^)^vi/^^^ = Q(g + G), 

ReA,B 

where A and B label the two sublattice, and, for R & A, we have taken by definition 
^i_R = ^i?, and '^2R = "^R+ax, being a the lattice spacing and x the unit vector in the 
x-direction. As soon as 7^ 0, nesting is not more important and indeed Q3 becomes 
massive. 
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Finally, because of GrGr and GaGa, also the conductance acquires other corrections 
with respect to standard localization theory. Indeed, these corrections add to give a vanishing 
/3-function for a^x-, as first indicated by Gade and Wegner [§. We have explained in the 
Introduction section (see also Appendix E) that this is a consequence of a simple abelian 
gauge symmetry generated by Eq.(|l|) at the particle-hole symmetry point E = Q. Similarly 
to the results for the density of states, this behavior of the /5-function is at odds with the 
standard theory. 



XIV. CONCLUSIONS 

In this work, we have derived the effective non-linear cr-model of a disordered electronic 
system on a generic bipartite lattice. This model, if the hopping matrix elements as well as 
the disorder only couple one sublattice with the other, shows an interesting behavior close 
to the band center, i.e. to the particle-hole symmetry point. Namely, the wave-functions 
are always delocalized at the band center, in any dimension. By a Renormalization Group 
(RG) analysis, in the framework of an e-expansion, e = — 2, we have found that the 
quantum corrections to the conductivity vanish if the chemical potential is exactly at the 
band center, thus implying a metallic behavior. In two dimensions, in particular, the Kubo 
conductivity flows to a fixed value by iterating the RG. On the contrary, we have found 
that the staggered density of states fluctuations, which are controlled by a new parameter 
in the non-linear a- mo del, are singular. This result is reminiscent of what it is found 
in equivalent one-dimensional models. In fact, models of disordered spinless fermions in 
one-dimension can be mapped, by a Jordan- Wigner transformation, onto disordered spin 
chains. In many cases, it is known that, in spite of the presence of disorder, these spin chain 
models display critical behavior, as shown in great detail by D. Fisher for random Heisenberg 
antiferromagnets and random transverse- field Ising chains |T^. Indeed, as pointed out by 



Fisher, the staggered spin fluctuations, in a random antiferromagnetic chain, also display 



critical behavior in the form of a power law decay, (— l)-^(5'(i?))(5'(0)) ~ i? ^, where the bar 
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indicates impurity average. Since the staggered spin-density corresponds to the staggered 
density of the spinless fermions, this resuh is consistent with the outcome of our analysis, 
which further suggests that a similar scenario generally holds in such models. Moreover, 
as in one-dimension [P,p|,p!7|, we find that the density of states is strongly modified by the 
disorder at the band center, and it actually diverges ind = 2. In reality, a random Heisenberg 
chain, away from the XX Z limit, maps onto a spinless fermion random hopping model 
in the presence of a random nearest-neighbor interaction. However, even in the presence 
of this additional interaction, the Hamiltonian has still the abelian gauge-like symmetry 
described in Appendix ^ which is at the origin of the delocalization of the band center 
state. This observation is also compatible with Fisher's result that the physical behavior 
does not qualitatively change upon moving away from the XXZ limit towards the isotropic 
XXX Heisenberg point. 

Many of the results which we have derived were already known. The existence of de- 
localized states at the band center of a two-sublattice model was argued already in 1979 
by Wegner The effective non-linear cr-model when the disorder breaks time- reversal 

invariance, as well as the RG equations, have earlier been derived by Gade |]^, although in 
a particular two-sublattice model. The extension to disordered systems with time-reversal 
symmetry was later on argued by Gade and Wegner P]. Finally, random flux models and 
disordered Dirac fermion models have recently been the subject of an intensive theoretical 
study [p!0|-[T2|, for their implications to a variety of different physical problems. 

In spite of that, our analysis has several novelties with respect to earlier studies. First 
of all, the two-sublattice model which we study is quite general. Secondly, the physical 
interpretation of the parameters which appear in the non-linear cr-model is quite transpar- 
ent. Thirdly, the explicit derivation of the RG equations with time-reversal invariance is 
presented. 
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APPENDIX A: SPECIFIC EXAMPLES 

As an example, we consider a tight binding model with only nearest neighbor hopping 
on a square and honeycomb lattice. 

In the case of square lattice, the enlarged unit cell is the x "\/2 one. The new reciprocal 
lattice vectors are Gi = 2n(l, —l)/a, G2 = 27r(l, l)/a, and the angle 9k of Eq.(||) is 



(Al) 



In the case of the honeycomb lattice, the unit cell contains already two lattice sites. The 
energy is given by 

1 + 2 cos (^-kxO^ cos i-^kya 



+ 



2sm ( -Kr^aj cos I —^kya 



and 



Ok = tan 



^ / 2 sin Q/cxo) cos {j^kyC 



1 + 2 cos (ffcj^a) cos {jY-ky 

The Brillouin zone is still honeycomb, with the y-axis one of its axes, and side equal to 
47r/(3V3a). 



APPENDIX B: WARD IDENTITY 



Let us consider a generic Hamiltonian in the two sublattice representation 



■2' 



(Bl) 



Ri,R2 
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where Hr^^r^ is a 2 x 2 Hermitean matrix. The current operator 

J{R) = E 

Ri,R2 

can be obtained by the continuity equation, leading to 

V ■J{R) ^ i^c^RHR^R^Cji^ - c^ji^HR^^RCji, 
Ri 



(B2) 



being V the discrete version of the differential operator. The long-wavelength expression for 
J{R) can be obtained by Fourier transformation, namely, through 

i^g • J(i?)e-'^-^ = i Y: {c^rHr,r,Cr^ - c^r,Hr,,rCr) e^^^'^, 

and expanding both sides in g, we get, for the linear term. 



^Y^Q-J{R)= E ^ • ^ {c^rHr,R,Cr^ - C^r^Hr.^rCr 
R RjRx 

J2 Q- {Ri - R) C^R^Hr^^rCr, 



R,Ri 



hence 



Ri 



Let us define the correlation functions 



JiR) = -i E (^1 - R) ci,HR,,RCR. 



(B3) 



(B4) 



where /i — 0,1,2,3, Jrj^^r^{R) — Srr^Srr^ are the density matrix elements, and J\ for 
i = 1, 2, 3, are the matrix element components of the current. By the continuity equation, 
we find that 

dtXo,i + djXj,i = ^ E -Sit - t^)Tr [G{Ru R\ ^2 - t)rR^R, {R'] 
Ri 

+6{t - h)Tr \G{R, t - h)J},^^j,{R') 



If we integrate both sides by 
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at = we find 



djXjAR^ R'; E)=tJ2Tr \GiR, R,- 



— Tr 



G{R,,R-E)J^^^^^{R') . (B5) 



Using once more the continuity equation ([B2|), we find 



d,d[x,AR^ R'; E) = -Tr [G{R, R!- E)Hr,,r + G{R!, R- E)Hn,R' 

+6Rn' E Tr [G{R, Ri; E)Hr,,r] + Tr [G{Ri, R; E)Hr^r,] . 
Ri 



By Fourier transform, 



Y: mxjAR, R'; E)e-^^^<^-^'^ = y: m,x,AR. R'; E)e-^i<^-^'^ 



RR' 



RR' 



Y[l- e-^^"(^^-«^')) Tr [G{R, R'; E)Hr,,r + G{R! , R; E)Hr,r, 



RR' 



At small q, the above expression is 



J2 (ii(ijXj,i{R^ R''^E) 



RR' 



- Y mj {R^ - Rd (Rj - R-) Tr [GiR, R'; E)Hr,,r + G{R\ R; E)Hr,r,- 



RR' 



leading to 



Y XjAR^ R'; e) = Y iR^ - R'i) {Rj - R'j) Tr [G{r, r!- e)Hr.,r] 



RR' 



RR' 



(B6) 



APPENDIX C: LONGITUDINAL MODES 



As discussed in Section |VII| , the expression of the Q-matrix which includes also the 
longitudinal modes is Qp{R) = Q{R) + S{R), where Q{R) and S{R) have been defined 
through (|7^) . The free action for these fields contains a local term, Eg . (|76D , and a non local 
one. Eg. (fTTp . The latter is 

1 



dRdR'T{R- R')Tr ArQ{R!)ArQ{R!)^ 



+ ^rS{R!)^rS{R!)^ + 2/\rQ{R!)^rS{R!) 



(CI) 
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where we have defined the operator 



Anf{R') = fiR) - f{R') = Y.^{R-R'T- V"/(i?' 



n=l 



Let us apply this operator to Q{R) and S{R), keeping all terms which contains at most two 
derivatives which act to the transverse matrices T. We obtain 

(C2) 
(C3) 
(C4) 



AnQ{R') ^ VQ{R') ■ [R - R'] 
AnSiR') ^ fiRyAnPiR')T{R') 



+ 



R-R 



+ 



Vf{R'y) P{R)T{R') + f{RyP{R) (yT{R'] 
V^f{R'y) P{R)T{R') + 2 (yf{Ry) P{R) (VT(i?' 



+ T{R!yP{R) N^T{R! 



R-R 



(C5) 



The term which is obtained by ( |CB| ) times its hermitean conjugate, together with the local 
piece (^6]) give the free action of the longitudinal modes 



(C6) 



The mixed terms, after defining D = TVT^ give rise to the coupling between transverse 
and longitudinal modes 

S[Q,P] = -\j dRdR'T{R-R') [R-R')Tr [d{R') (^P\R')P{R) - P{Ry P{R'))] (C7) 



dRdR'r{R- R') (^R-R'Y 



Tr 



aiD{R')aiP{R)D{R')P\R') + aiD{R')aiP{R')D{R')P\R) 



D{R')D{R!) (P\R)P{R!) + P\R')P{R) 



^ / dRdR'T{R - R!) [R - R') ■ Tr VQ{R')ARP{R'y + H.c. 



(C8) 
(C9) 
(CIO) 



The last term, ( |C1CI|) , gives rise to higher gradient contributions, hence can be neglected. 



1. Longitudinal propagators 



Before averaging over the longitudinal modes, we have to evaluate the longitudinal prop- 
agators. The matrix 
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P = {PofiSo + Po,3S3) CTq + i (Ps.lSl + -P3,2S2) CTg, 

has to satisfy cP*c* = P, and, in addition, all Pa = P]. For a = (0, 0), (0, 3), (3, 1), by 
writing 

we find that 

pf,p.e7^e, pr = (pry, p.=-(p.)* 

For q; = (3, 2), by writing 
we must impose 

-^3,2 ) -^3,2 fc /v-e, -r3_2 — ~ 1^-^3,2 J ) -^3,2 — [^3,2 J ■ 

If P^^ is a symmetric real matrix, its propagator is 

= f (<^a.<^6c + , (cii) 



while if it is antisymmetric 

(^i:L<L) = - ^ (5ad<^6c - SM , (C12) 



where Gk — cuk/SV, and {a,b,c,d) are replica indices. By means of these propagators, we 
readily find that, if M = M^^^ + iMj • r, where M*^*) are matrices in the replica space, then, 
for a — (0, 0), (0, 3), (3, 1), the following results hold 

{PaMPa) = -GcM^c* + 2GTr (m^^^) = -GcM^c^ + GTr (M) , (C13) 

while, for a — (3, 2), 

(P3,2MP3,2) = GcM^c* + 2GTr (m(°)) - GcM^c' + GTr (M) . (C14) 
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More generally, 

{PMsiajP^) = -GcM^c^ {siaj + SsSjCTjSa + crsSiSjO-jSicrs - cr3S2Si(7j 82(73) 

+ GTr (M) {si(jj + ssSiajSs + asSiSiajSias + 0-3525^0-^520-3) , (C15) 
{PMsiajP) = -GcM^c^ {siaj + 535^0-^53 - o-3SiSiO-jSiO-3 + 0-3325^0-^520-3) 

+ GTr (M) {siaj + 535^0-^53 - a^SiSiajSias - a3S2SiajS2cr3) . (C16) 

For j = 0, 3 the above expression simplifies to 

{PMsiajP^) = -2Gc {Msiajf c* + 2GajTr (Ms^sq) , (C17) 
(PMsiajP) = -2GS3C {MsiGjf c^ss + 2GajTr (MsiSa) , (C18) 

while for j = 1,2 

(PMsiajP^) = -2GS3C {MsiajY c's^ + 2GajTr (MsiS^) , (C19) 
{PMsiajP) = -2Gc {MsiajY c* + 2GajTr (MsiSo) . (C20) 



2. Averaging S[Q,P] 

We have now all what it is needed to proceed in the averaging over P. Here we just sketch 
the calculation, which is quite involved and requires the matrix properties of D which are 
determined in Appendix We just remark that ( |C8| ) and ( |C^ ) do not reproduce the 
correct stiffness term. Indeed, it is (|C7|), which contributes at second order, which cancels 
the additional terms and allows to express everything in terms of the matrix Q. 

By means of the previously calculated propagators of the longitudinal modes, we find 
that 

{S[Q,P])p = ^ I dRTr[VQiRyVQiR)]+^[Tr[QiR)^VQ{R)a3)]\ (C21) 
where 

Y = j dRT{R)T{R)-^R^. (C22) 
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We notice that the first term is a contribution to the Kubo conductivity and the second 
to the staggered density of states fluctuations of the diagrams where the current vertices are 
those proportional to the random hopping. 



3. Additional terms 

The last class of corrections which generate new operators is obtained by expanding each 
Green's function in (|50|) at first order, Gp = Gq — iGqPGq, leading to the term 

^{Tr{GPGUGPGU))p. (C23) 

For the sake of clarity, we will analyse this term only in the case of a cubic lattice, where the 
derivation is more straightforward. We will postpone a discussion about the general case at 
the end of the section. 

In the cubic lattice, according to Eg. , the electron-Q coupling can be brought to a 
local one also in the diagonal basis. In this basis, Q = Qq + iQicri, 

PiR) = (Pooso + P03S3) c^o + ^ (Pusi + P12S2) (C24) 
the unitary matrix 



T{R) = exp 



Wo{R) , W,{R) 



(C25) 



2 " 2 

where Wi has the same form as W3 defined by Eg. (|3^) , and T = a^Ta^. Moreover, since 
Hrr' = eR-Ria^ in the diagonal basis, being eR-Ri the Fourier transfom of e^, the current 
operator which appears in the definition of Urri, Eg.(^), has only matrix elements Jk = 
Vefc 0-3. 

Once averaged over P, (|C23|) gives, among other terms which correct the Kubo conduc- 
tivity, a new term [see Eg.(Greenfunctiond)] 

a, 6=1, 2 h=± pkq 
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where the first piece derives from Pq and the second from Pi. The structure in the en- 
ergy/sublattice indices can be shortly represented by 

Aj (g) Ai = ^ (a-2 ® 0-2 + S30-1 (g) saCTi + 0-3 (g) (T3 + S3 ® S3) , 

j=l,...,4 ^ 

SO that the above term can be written, at small q, as 

I E J2^p-kTriUggpAigp)TriU-ggkAigk) (C26) 

j=l,...,4 pkq 

We remind that 

Gp^qGp ^gp-Dq ■ Jpgpj 

where Dg is the Fourier transform of T{R)'VT{R)~^ . We notice that only the diagonal 
component in energy of D enters. Moreover, for the diagonal matrices Sj = Sq, S3, since 
D = -Dqcto + -DiO"i, it derives that, for i = 0, 1, the following equality holds TrDsjai = 
—Tr {vV/sjOi/'l^. Therefore, just Wiai contributes. By means of (pTl) , for any A's we have 

-iTr ipi^qai ■ Jpg{p)Aig{p)) = +Vep • Tr (-Di^<ySi^(p)sicr2Ai^(p)) 

= Vep ■ Tr [Di^g [Gla3a2Aia3 - Gls3a2AiS3 + G'3G'oS3[a3, a2Ai])] . (C27) 

Only for Ai = (72/4 the trace over the a's is finite, leading to 

2 (gI - Gl) Vep ■ Tr J . (C28) 



In conclusion, going back to the original sublattice representation, and defining 

n(°) = E Vefc ■ Vep Up_,Tr (a3G+a3G+) Tr (a3G+a3G+) , (C: 

pk 

we obtain the following explicit expression of ( |C23|) [notice that Q = QoCo + iQiO'i 
Qo<^o + iQsO's in the sublattice basis] 
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n(°) J dRTr [q\R)VQ{R)(T3] ■ Tr [Q\R)VQ{R)a3] . (C30) 



8 ■ 32S4 

This term, which derives from the regular current vertices, has the same form as the second 
term in ( |C21D , which, on the contrary, is due to the random current vertices. Therefore, the 



coupling constant 11 which appears in the final expression (|80D is the sum of both terms. 
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The symmetry of the operator ( |U30|) , which, due to the trace, involves the Nambu, 
energy and sublattice components tq, sq and o"3, respectively, suggests that the prefactor 
n represents the fluctuations of the staggered density of states, as discussed in section 
Vllllj| . In the case of a generic bipartite lattice, we do expect an analogous term to appear, 
because the staggered-density of states fluctuations still acquire singular contributions, and 
the operator is not forbidden by the symmetry properties of the Q-matrix. The reason 
why we decided to show only the case of cubic lattices is that the distinction between the 
longitudinal from the transverse modes is a bit ambiguous at large momenta, where the 
both are in a sense massive. This is not a problem for cubic lattices, where one can show 
that the large momentum components of the transverse modes do not contribute, hence 
( P30|) exhausts the whole contribution. On the contrary, in other cases, we do have to keep 
into account the contribution of the small- wavelength transverse modes to recover the full 
expression. This makes the calculations more involved than in the case of cubic lattices. 



APPENDIX D: EXPLICIT DERIVATION OF THE RG EQUATIONS 

In this Appendix we outline the derivation of the RG equations. Before that, it is 
convenient to list some useful results. 



1. Averages 

Besides the propagators, we will also need the explicit expression of particular averages 
which enter in the derivation of the RG equations. The following results hold: 

{BabPbcBl) = W{k)5adTr (Po) , (Dl) 
{BabPbcBcd) = -2Dik)Pl, (D2) 

where P = Pq + iP ■ r is a quaternion real matrix. 
In addition. 
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{■^abPbc-^cd) — -2D{k)Pl, + AD{K)5aaTr{Po) 



-m- 



n 



' + rim" 
{CabPbcCcd) = -'iD{k)Pl + 4L>(K)5„<irr (Pq) 



D{k) 



n 



(J^a; + rim" 



n 



a^^ + Hm 



For instance, if P = 7, then 

(PS^) = AD{k)mi 
{BB) = -2D{k)i 



{CC) = 



-2D{k) + AD{k)m - D{k) 
-2D{k) + 4D{k)m - D{k) 



n 



a XX + rim 

n 



a XX + Ilm 



(AC) = -D{k) 



n 



(J^a; + Ilm 



(D3) 

(D4) 
(D5) 



(D6) 
(D7) 
(D8) 

(D9) 

(DIO) 



2. RG equations 



First of all, we need to know the quaternion structure of the matrix D = TVT^ . Since 



D = -D\ cD^c^ = -aiDai, 



we can write D = P'oCo + -Dscxa, where in the ±-space 

Do 



^ Ao Po^ 



V 



and 



P.. 



Bl Co 

A3 ^3 
Bi C3 



(Dll) 



Wc can write each of the above matrices in quaternion form, P 
p(o)^ P e TZe, but, in addition, we must impose that 



(D12) 



P^^^tq + iP ■ t, where 



50 



By making use of the above properties, and by means of the Eqs.( P2|) , (p3|) , 
after defining 



41) and 



A/s (27r) 



we get the following results 

{WfDoWf)f = LATDo 
+2Lb 



/ 


Bo) 




%o 


\ 




-2La 




V 


-Bl j 




^ Co 





(D13) 



where F = 11/ {a^^ + mU), and 

{WfD,Wf)f = LArD3 









f 




+2LAi 




- 2LBi 













-ALITr iA 



l(0) 



(0) 



(D14) 



For further convenience, when useful, we have labelled Lb the propagator of i?/, and La the 
ones of and Cf. The next useful result is 



{WfWf)f = {-ALbtti - ALatti + 2La + LaT) L 
Through ( P15 ) and (P14 ) we therefore get 



(D15) 



JdR- 2{Tr [DaiQspWfDWfQ,pa,])f 
+2{Tr [Dcj,QspDQspWfai])f 

271(7 XX 
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j dRA{LB - La + 2LAm + 2LBm) Tr 



I Bo^ 



\-Bl J 



+ 4 (2LA'm + 2LBm) 



^^3 ^ 

iC^ 
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(D16) 



(D17) 



-2L[Tr {Da^)f 

- 4 {2LAm + 2LBm + 2La) 



Co 
/ 



- 4 {2LAm + 2LBm + La + L^) 





z5l 



(D18) 
(D19) 

(D20) 



The calculation of (-Ff )/ is more involved, since one needs the average of four VF/'s. For 
sake of lengthy, we just quote the final result that such a term cancels ( pi9| ) and ( p2U| ). We 
next notice that 

1 



— Tr 
E2 



VQ^VQ = 2Tr [Ds^aiDs^ai - DD] 

2 



-4Tr 



/ 


Bo\ 


2 




^ lA, ^ 










-4Tr 




V 


-Bl j 






^ ^C3^ 



and that 



Therefore, ioi La = Lb, 



Tr 



gtygasl = 2S2rr {Da;) . 



= ^4Lm'^^^ I dRTr 



32E2 



1 ^ 27ro"^.T 



(D21) 



VQ^VQ 

2 32S4 

In the cases in which the A and C modes are gaped {La = 0, and no D^), we obtain the 
standard result 



{F,-lF^)f = -{2Lm + L)^ I dRTr 



(D22) 



APPENDIX E: GADE AND WEGNER'S PROOF OF THE VANISHING 

/3-FUNCTION 



The equations from ( |5^ ) to (^) imply that W; is not a traceless matrix. Indeed, we can 
alternatively write W as 
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W = W' + -^Tr (Ws) a3 = W' + i(pa3, (El) 
4m 

with W = Woao + W^a^, now being a traceless matrix. Since cxs commutes with W, 
this means that 

T{Ry = e^(^) = e^<^(«)'^«e^'(^) = e^'^^^'^'^W (R) , (E2) 

which also defines the matrix field V{R). By means of this parametrization, the non-linear 
cr-model (j^) can also be written as 

S[T] = S[V, (j)] = / dR Tr ^V{R)-^ ■ VV{R) 

32 J 

7T17T f 

+ -Y (^^^ + j dR V(f){R) ■ V(f){R). (E3) 

Therefore the action of V is distinct from that of (p, and the latter, being a phase, is gaussian. 
This implies that the combination axx + is not renormalized and scales with its bare 
dimension e, for any number of replicas. In turns, it means that, in the zero replica limit, 
it is which is not renormalized! This is completely equivalent to the nice proof given by 
Gade and Wegner that the quantum corrections to the /3-function of the conductance of 
a U(A^)/SO(A^) model vanish at all orders in the ^ limit. 

The other important result concerns the renormalization of an operator 



Within RG, 



^iq^a^ ^ -g2 In s i —\ Q^l't'-i . (E4) 

2 V c + mt m J 

The second term, which is singular in the m — limit, has to be canceled by the one- loop 
renormalization of V^. Gade and Wegner showed that this cancellation holds for any q. 
Furthermore, they argued that, apart from the one-loop correction, the renormalization of 
does not contain any other singular term in the m — > limit. This, as they pointed out, 
has very important consequences. In 2d, t does not scale, while c goes to zero. Therefore 
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the term which dominates the renormahzation of T^'' for m = is just the first term in the 
right hand side of ( |E4| ) . This argument imphes that the one-loop correction, which we have 
derived for the density of states {q = 1 case), is sufficient to identify the correct asymptotic 
behavior. 

To conclude, let us discuss more in detail the origin of this gaussian field 0. In the 
Grassmann variable path-integral representation, the action for the particle-hole symmetric 
model at E = u = 

S = ^"^RHnR'^i/R', 

RR' 

posseses a simple abelian gauge-like symmetry 

^ ^ e^'^'^^^, (E5) 

because {as, Hrri} = 0. It is just this symmetry which causes the appearance of the gaussian 
part of the non-linear cr-model. Notice that this symmetry implies a particle-hole symmetric 
Hamiltonian, which is invariant under the transformation ^ — > cj ^ but C2 ^ ^ —cl ^. In 
fact, {(Ts, Hrri} = also means that {ticts, Hrri} = 0, being Hrr/ oc tq. The interesting fact 
is that ( |E5| ) is not a symmetry of the fermion operators. Indeed, under this transformation, 
cn — >• e^^'^'-^CR, but cr e^'^'^^Cij, and not cr — >• e~'''^''^CR as we would expect if cr and cr 
had to be identified with the operators Cj^ and c|j. Finally, we notice that if, besides (T3, the 
Hamiltonian commutes with another Pauli matrix (as it can be the case for specifically built 
particle- hole symmetric models), the above gauge symmetry would be non abelian, hence 
spoiling all peculiar properties which we have shown to occur. Indeed, in the last case, 
the system can be mapped into a standard localization model with an additional sublattice 
index. 
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